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Abstract 

CN ■ Solution of the Schrodinger's equation in the zero order WKB approximation 

. is analyzed. We observe and investigate several remarkable features of the WKBq 

^SJ I method. Solution in the whole region is built with the help of simple connection 

■ formulas we derive from basic requirements of continuity and finiteness for the 

^ . wave function in quantum mechanics. We show that, for conservative quantum 

systems, not only total energy, but also momentum is the constant of motion. We 
derive the quantization conditions for two and more turning point problems. Exact 
energy eigenvalues for solvable and some "insoluble" potentials are obtained. The 
eigenfunctions have the form of a standing wave, AnCos{knX + 6n), and are the 
asymptote of the exact solution. 
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^ 1 Introduction 

O I The wave equation in quantum mechanics can be derived with the help of the Bohr's cor- 

respondence principle. This fundamental principle has been used at the stage of creation 
of quantum theory. It is used to establish correspondence between classical functions 
■ and operators of quantum mechanics, and to derive the apparent form of the operators 

^ . [0) 0- Moreover, the correspondence principle points out the way to a simplest asymptotic 

Q-i| solution of the Schrodinger's equation. 

The correspondence principle states that the laws of quantum physics must be so 
chosen that in the classical limit, where many quanta are involved, the quantum laws 
lead to the classical equations as an average. In this way, in Ref. this principle has 
been used to derive the nonrelativistic semiclassical wave equation appropriate in the 
quasiclassical region and, in Ref. 0] , the corresponding relativistic wave equation. It was 
shown that application of the standard WKB method to the semiclassical wave equation 
allows to solve the so called "insoluble" problems [^, H, |[ and results in the existence 
quantum fluctuations of the angular momentum, which contribute to the energy of the 
ground state 0. 

There is a mathematical realization of the correspondence principle known as the qua- 
siclassical approximation which is widely used mainly as the Wentzel-Kramers-Brillouin 
(WKB) method [§, [ill applicable in the case when the de Broglie wavelength, 

X = h/p {h = 2TTh), is changing slowly. This method is usually used as a tool to obtain 
the approximate solution of the one-dimensional Schrodinger's equation in the quasiclas- 
sical region at large values of quantum numbers. 

In this work, we analyze the classical limit of the Schrodinger's wave mechanics with 
the help of the correspondence principle. We investigate solution of the wave equation 
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in the zero-order WKB approximation (WKBq) and observe several remarkable features 
and important advantages of the zeroth approximation with reference to the commonly 
used first-order WKB approximation (WKBi). These advantages are: 1) the WKBq 
wave function (w/), iPq{x), has no divergencies at the turning points; 2) the WKBq ap- 
proximation allows a considerably simpler (than the standard approach based on linear 
approximation) derivation of the connection formulas and quantization condition; 3) the 
WKBq method allows to solve the wave equation for many-turning-point problems; 4) for 
stationary states, not only energy, but also momentum is the constant of motion; and, 
finally, 5) the most important advantage is that the WKBq approach reproduces the exact 
energy spectra for all known solvable potentials and "insoluble" potentials as well. 

Solution of the wave equation in the whole region is built with the help of connection 
formulas we derive in this work, which give the main term of the asymptotic series. To 
demonstrate efficiency of the WKBq approach, we calculate exact energy eigenvalues for 
solvable and some "insoluble" spherically symmetric potentials. For discrete spectrum, 
the WKBq eigenf unctions are written in terms of elementary functions in the form of a 
standing wave, y4„cos(A;„x + and give the asymptote of the exact solution. 



The zero order WKB solution 



It is well known that the exact eigenvalues can be defined with the help of the asymptotic 
solution, i.e. the exact solution and its asymptote correspond to the same exact eigenval- 
ues of the problem under consideration. The asymptotic solution in quantum mechanics 
can be obtained by the WKB method. This means that the quasiclassical method can be 
considered as a tool to reproduce the exact energy spectrum. This is true for some poten- 
tials [jl2|, |13[ 0, and has proven for many others with the help of specially developed 
techniques, improvements, or modifications of the quasiclassical method on the real axis 
and in the complex plane 
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The main problem concerning the applicability of the WKB approximation is the 
problem of exactness of the method. Proofs of varying degrees of rigor have been ad- 
vanced that demonstrate the exactness of the standard WKB quantization condition 

The question of exactness of the WKB 
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approximation is usually reduced to the estimation of the high-order correction terms. 
The earliest development of the WKB method for obtaining the high-order corrections 
has been considered in Ref. 
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Then, in the radial generalization of Dunham's 
one-dimensional WKB quantization condition was derived with the help of the Langer 
transformation UTS). It was shown that the second- and third-order integrals identically 
vanish for the hydrogen atom and the three-dimensional harmonic oscillator |13[. As we 
show below, the first order term also does not affect the exactness of the quasiclassical 
quantization condition. 

Consider the static Schrodinger's equation with the arbitrary potential V{x) in one 
dimension, 




dx 



V{x) 



ip{x) = Eip{x). 
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Let us write the equation (|T]) in the form, 



-ih— -^(x) = - Uix) ^pix) 
ax L -I 



(2) 



where = 2mE and U{x) = 2mV{x). Equations (|I]) and (^) are equivalent, but the 
physical sense of their solutions is different. The Schrodinger's equation (|l]) is the energy 
eigenvalue problem, and equation @ is the eigenvalue problem for the squared momentum 
P^. Equation (^) has two solutions for each eigenstate (as it should be for the second- 
order differential equation) with the corresponding eigenvalues P^''^^ = ^\/2mE^^ i.e. 
eigenvalues for these two problems are connected via the equality P^ = 2mEn for a free 
particle. 

In the conventional WKB method, solution of Eq. (0) is searched for in the form 



x) = Aexp 



-S{t,x) 



(3) 



where A is the arbitrary constant and S(t,x) is the complex function. Substituting (|) 
into Eq. (Q) one finds that S{t, x) satisfies the equation 



^ dS\ fi'^ ( p2 
dx I dx \dx I 



U{x). 



(4) 



The Hamilton- Jacobi equation in quantum mechanics (^) is equivalent to the wave equa- 
tion (H). In the limit h —>■ Eq. reduces to the classical Hamilton- Jacobi equation. 

In the WKB approximation one expands S(t,x) in power series of h, S = So+ hSi+ 
h'^S2 + . . . and usually keeps two first terms of the expansion, i.e. 5*0(^,0;) and Si{t,x) 
(the WKBi approximation). This gives, after simple calculations, (dSo/dx)^ = P^ — U{x), 
and Si{t,x) = iln ^J\p{x)\, where p{x) = dSo/dx. For conservative systems, the classical 
action So(t,x) can be written as So(t,x) = —Et + W{x), where W{x) is the Hamilton's 
function (reduced classical action). Hence, one gets 



'dW^ 
dx 



P^ - U{x). 



(5) 



Here dW/dx = p{x) = ii/P^ — U{x) is the generalized momentum and the WKBi wf is 



(1) / 

WKB\ 



where 



\pix)\ 

[x) is dimensionless phase variable. 



(6) 



-W{x) 



^P"^ - U{x)dx . 



(7) 



1 

h" ~h 

There are several disadvantages of the WKBi approximation. The most essential of 
them is that the WKBi wf (^ diverges at the classical turning points given by p{x) = 
0, where the condition for its applicability breaks down. Although this divergence is 
understandable in the classical limit, since a classical particle has zero speed at the turning 
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points, it is certainly not present in a quantum mechanical treatment. Since "ipwKBi^) 
singular, tricky matching the wf ipi^j^ ^{x) and its first derivative at the classical turning 
points have been developed [1^ that yields the famous connection formulas and well 
known WKB quantization condition. The exact mathematical conditions under which 
the connection formulas can be proved rigorously to apply are fairly complex. Other 
difficulties arising in the WKBi approximation have been considered elsewhere 0. 

The connection formulas and quantization condition can be obtained much simpler 



in the WKBq approximation I^O], |21|], i.e. in the classical limit ^ ^ 0. In this case 
S{t, x) ~ So{t, x) and the general solution of the wave equation for the potential V{x) 
can be written in terms of the dimensionless phase variable 4>{x), 

ij^^j^ix) ^ ^o[0(x)] = Cie^^(^) + C2e-^<^(^') , (8) 

where Ci and C2 are taken constant in the classical limit. 

The WKBq solution (|^) has several important advantages with reference to the com- 
monly used WKBi wf(^. The most important is that the solution has no divergence 
at turning points. In the phase space, the function (H) has the form of superposition of 
two plane waves; it has also such a form in the configuration space (see below). Note, 
that the WKBq approximation corresponds to the main term of the asymptotic series in 
the theory of the second-order differential equations. 



3 Connection formulas. Quantization 

Basic requirements for the wf in quantum mechanics are continuity and finiteness in the 



whole region. The famous connection formulas [0] have been derived to connect the 
WKB solution in the classical region with the one in non-classical region. These formulas 
are based on a linear approximation to the potential in the regions around the classical 
turning points where the WKBi w/ diverges; the exact solution involves Bessel's functions 
of order 1/3. Even though the connection formulas themselves turn out to be simple 
enough, their derivation is quite tedious. 

One of the most important problems to which the connection formulas apply is that 
of penetration of a potential barrier. Interaction with the potential barrier can be of two 
types - refiection and transmission. There are several problems in which the conventional 
WKBi formulas break down [||]. One kind of problems occurs near the top of a barrier, 
where the slope of the potential is small. As a result, the straight-line approximation to 
the potential breaks down, and the connection formulas must be altered. The second type 
of problems arises when the potential changes too rapidly, for example, in the case of a 
square well potential. The connection formulas also break down for some other type of 
potentials [Q]. 

Consider the WKBq wf (|^) for the two-turning-poing (2TP) problem. To build the 
physical solution in the whole region we need to connect the oscillating solution in clas- 
sically allowed region with the exponentially decaying solutions in classically inaccessible 
regions. For the 2TP problem, the whole interval {—00, 00) is divided by the turning 
points Xi and X2 into three regions, {—oo,Xi] (region I), [xi,X2] (region II), and [0:2,00) 
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(region III). In the classically allowed region II, right from the turning point Xi, the real 
part of solution can be written in the form 



V^o''(0) =C^icos(0-0i + (5i), (9) 
and left from the turning point X2, 

4'i<l^)=C2COs{<j)-<f)2 + S2), (10) 

where 0i = 4>{xi), 02 = 0(2^2), and 61^2, are the phase shifts at the turning points xi and 
X2, respectively. In the classically inaccessible regions I and III, to guaranty finiteness of 
the wf at X ^ ±00, we choose the exponentially decaying solutions, i.e., 

4{<P) = Ae'f'-'^\ (11) 

left from the turning point Xi, and 

^o'^^(0) = Be-^+<^^ (12) 

right from the turning point X2- 

To construct solution in the whole region we need to connect functions ipQ^cf)), 
and ipQ^\(j)) at the turning points Xi and X2- Matching the functions ipQ{<f)) and ipl^cj)) 
and their first derivatives at the turning point x = Xi gives 



J Ci cos 5i = A, 
1 Ci sin 5i = —A, 

and for ipQ^{4>) and ipo^^i^') at the turning point x = X2 we have 

f C2 cos(52 = B, 
1 C2 sin 52 = B, 



which yields 



and 



tan^i = — 1, 
tan52 = 1, 



Ci = V2A, 
C2 = V2B. 



(13) 



(14) 



(15) 



(16) 



The connection formulas (|T^) and (|T6D supply the continuous transition of the oscil- 
lating solutions (I) and ( |T0|) into exponential solutions (pTD and ([T2|) at the turning points 
xi and X2. 

The function ipl^cf)) [given by Eqs. (^ and (0)] should match itself at each point of 
the interval [xi,X2]. Putting x = X2 we have, from Eqs. (|^) and (pUj), 

CiCOs(02 - 01 + 5i) = C2C0S52- (17) 

This equation is valid if 
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02 - 01 + 5i = ^2 + 7rn, n = 0,1,2,... (18) 

and C2 = (— l)"Ci. Equation (^) is the condition of the existence of continuous finite 
solution in the whole region. This condition being, at the same time, the quantization 
condition. Taking into account the notation (|^, we have, from Eq. (0), 

^P"^ - U{x)dx = Tin[n + ^. (19) 

Condition ([T9|) solves the 2TP eigenvalue problem given by Eq. (^. We see that the 
WKBo approximation results in the quantization of the reduced classical action, W{x), 
and coincides with the well known WKB quantizaton condition. It is important to under- 
line that the first order term does not affect the exactness of the quantization condition. 
In our case, the condition (^^ is obtained by product from general requirements of conti- 
nuity and finiteness to the wave function in the whole region, i.e. from the requirements of 
a smooth transition of the oscillating solution given by Eqs. (|^) and (p^Of ) to the exponen- 
tially decaying solutions (pJ]) and (|12D in the classically inaccessible regions. Quantization 
condition ([ToD reproduces the exact eigenvalues for all known 2TP problems (see Refs. 
i I, 0)- 

The derivation given above is evidently much simpler than the usual text book ap- 
proach for deriving connection formulas and quantization condition. In the general case 



of three or more turning points, phase shifts 5i and 82 will be different . 

Combining the above results, we can write the finite continuous WKBq solution in the 
whole region (^2 = — = vr/4), 

V^o[0(a;)] =c| cos[0(x)-0i-f], xi < x < X2, (20) 
i ^^e-'^(^)+'^2^ X > X2. 

Oscillating part of solution (^) corresponds to the main term of the asymptotic series in 
theory of the second-order differential equations. In quantum mechanics, the oscillating 
part of Eq. (pO|) gives the asymptote of the exact solution of the Schrodinger's equation 
(see below). 



4 Momentum conservation for discrete spectrum 

It is well known that, for discrete spectrum, the asymptote of Eq. (0) has the form 
ipni^) = AnCos{knX + 6n), whcrc A„, /c„ = Pn/h, and 6n are constant values; here A;„ is 
the wave number and P„ is the momentum eigenvalue. The asymptote has the form of a 
standing wave and can be written as a superposition of two plane waves, exp{iknx) and 
exp{—iknx) , each of which describes free motion. Therefore, for the discrete spectrum, 
the asymptote of Eq. (0) describes free motion of a particle-wave in the enclosure and 
the eigenmomentum P„ = hkn is the constant of motion. 

The quasiclassical approximation is the asymptotic method, which takes advantage of 
the fact that the wave length is changing slowly, by assuming that the wf is not changed 
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much from the form it would take if the potential V was constant. The WKBg solution (D 
is justified in the classical limit h —>■ 0, which corresponds to the short-wave asymptote 
of the Schrodinger's equation (0). Our basic requirements to the solution (H) are: 1) 
Ci and C2 are constant values and 2) the generalized momentum p{x) = dW/dx is the 
adiabatically slowly changing function, i.e., 

dW 

— — ~ const. (21) 
dx 

These are the same assumptions which have been originally used by Schrodinger to derive 
the wave equation from the optical-mechanical analogy (see, for instance, Ref. [jl|). 

In this analogy, the most important aspect is the relation to the limiting case of 
wave optics, i.e. geometrical optics. The limiting case of geometrical optics is analogous 
to the classical limit of quantum mechanics: the amplitude of the (wf) is a constant 
and momentum is the adiabatically slowly changing function. Note, the constraint (pl|) 
supplies the Hermiticity of the operator = [—ih{d/dx)]'^ in Eq. (g) 

Consider the WKBq eigenfunction (pOD in the classically allowed region given by p{x) > 
0. Show that the oscillating part of the wf (|3) is in agreement with the asymptote of 
the corresponding exact solution of Eq. (^ and can be written in the form of a standing 
wave. In equations (|l|) and @, the total energy E and the momentum P are connected 
by means of the equality = 2mE for a free particle. Integrating (]2T|), we obtain, for 
the action W{x), 

W{x) ~ Px + const, (22) 



where P = y2rnE is the total momentum. Solution of Eq. (|l|) [or (|^)] results in the energy 
quantization, i.e. E = E^. Therefore, W{x) ~ P^x + const, where P„ = hkn = \p2/mE^. 
Underline, this form of the action is valid only for allowed motions in quantum mechanics, 
i.e. for the energies E = En, which can be found with the help of the quantization 
condition (p!9|). 

Thus, the oscillating part of the wf (EOf) takes the form 



4\x) = Cn COS ( ip„x +\n]. (23) 



The normalization coefficient, 



Cn 



^ T:h{n + \) + h' 
is calculated from the normalization condition 



2P„ 



(24) 



\i)Q{x)y dx = I. (25) 

In Eq. (P^D, we have took into account the fact that, in the stationary states, the phase- 
space integral (|^ at the TP Xi and X2 is 0i = — 7r(?2 + l/2)/2 and 02 = 7i"(^ + l/2)/2, 
respectively, so that 02 — 0i = 7r(n + 1/2) [|, The phase integrals 0i and 02 depend 



on quantum number and do not depend on the form of the potential. This form of 0i and 
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(f)2 guaranties that the eigenf unctions are necessarily either symmetrical (n = 0,2,4, . . .) 
or antisymmetrical [n = 1, 3, 5, . . .). 

Solution (|23|) has the form of a standing wave, which can be written as the superposi- 
tion of two plane waves. The fact that, for stationary states, the oscillating part of solution 
( p3D can be written in such form means that (for conservative quantum-mechanical sys- 
tems) the particle momentum is the constant of motion and the corresponding coordinate 
is cyclical. Therefore, solution ( ]23|) describes free motion of a particle-wave in the en- 
closure, where the enclosure is the interaction potential. This means that interaction of 
the particle-wave with the potential reduces to reflection of the wave by the walls of the 
potential. 

As other integrals of motion in quantum mechanics, the total momentum P can take 
only discrete values, P = Pn, corresponding to the discrete values of the classical action 
W{x), which can be found from the quantization condition (|T^. Note that, for a particle 
moving in the Coulomb potential, the energy eigenvalues can be written in the form of 
the total energy for a free particle, i.e. En = Pn/2m, where P„ = iam{nh)~^ is the 
momentum eigenvalue and n is the principal quantum number. It is easy to show that 
the energy eigenvalues for other solvable potentials also can be written in such form. 

There are the following physical grounds to the above discussions. According to classi- 
cal theory, accelerated electrons radiate energy at a rate equal to |e^c"^x^. But electrons 
in atomic orbits are always accelerated (?) and should, therefore, lose energy continuously 
until they fall into the nucleus. Actually, it is known that this never happens. To resolve 
this problem, Bohr postulated that 1) the electrons are in stationary states, i.e. do not 
radiate, despite their acceleration and 2) electrons can take discontinuous transitions from 
one allowed orbit to another. However, the postulate is not an answer to the question: 
why the accelerated electrons do not radiate? 

Acceleration means change in momentum. In quantum mechanics, the change in 
momentum (which is associated with the de Broglie wave X = h/p, h = 2tt%) appears 
as radiation of frequency u = {E — E')/h that means change of the stationary state. 
As we have shown above, in the stationary states, particles are not accelerated, because 
the momentum is the constant of motion. This means that, in the stationary states, the 
electrons move as free particles, particles-waves in enclosures. 

The well known quasiclassical condition. 



p 



dp 



dx 



« 1 , (26) 



defines the region of applicability of the WKBi approximation. Conventional treatment 
of this condition implies that the momentum p{x) is large enough, i.e. quantum number n 
take large values. At the same time, the WKB method yields the exact eigenvalues for all 
values of n, including n = (see Ref. 0, ^, |T^). This means that the condition ( PBj ) can 



be treated differently: it supplies the hermiticity of the squared momentum operator, p"^, 
that implies dp/dx ~ 0, and does not imply that the momentum p{x) takes large values 
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5 The multi-turning-point problems 



The WKB method is usually used to solve one- dimensional two turning point problems. 
Within the framework of the WKB method the solvable potentials mean those potentials 
for which the eigenvalue problem has two turning points. However the WKBq solution 
(ID is general for all types of problems. The zeroth WKB approximation allows a consid- 
erably simpler derivation of the connection formulas and the corresponding quantization 
condition not only for two but, also, for many-turning-point problems (yuTP, /i > 2). This 
is the class of the so-called "insoluble" problems, which cannot be solved by standard 
methods. 

The quasiclassical formulas can be written both on the real axis and in the complex 
plane. Most general form of the WKB solution and quantization condition can be written 
in the complex plane. In the complex plane, the 2TP problem has one cut between turning 
points Xi and X2, and the phase-space integral (|r^) can be written as the contour integral 
about the cut. The fiTP problems contain (in general case) bound state regions and the 
potential barriers, i.e. several cuts. The corresponding contour C should enclose all cuts. 
The WKBq wf in the whole region can be built similarly to the 2TP one considered above. 

Let the problem has u finite cuts, i.e. ^ = 2v turning points, and the potential is 
infinite between cutJil. Then the integral about the contour C can be written as a sum of 
contour integrals about each of the cut, 

i _ u{x)dx = Y.j ^P"^ - U{x)dx, (27) 
where each term of the sum is the 2TP problem, 

£ y^p2 _ u{x)dx = 2Txh + ^ j . (28) 



Therefore, the /xTP quantization condition is |T9 



J C \ ^ J i 



(29) 



The condition (|29| ) is in agreement with the Maslov's theory |2^. This means that 
the right-hand side of the equation (P7| ) can be written as 

2r^hj2{n + ]^=2nh[N+^^. (30) 

where N = Yli=i is the total number of zeroes of the wf on the u cuts. In this interpre- 
tation, the number of turning points n is the Maslov's index, i.e. number of reflections of 
the wf on the walls of the potential. Thus the condition (|29|) takes the form 



£ _ u{z)dz = 2Txh [n + (31) 



^In case where the potential is finite in the whole region, the quantization condition will be more 
complicate p2[. 
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6 Multi-dimensional problems 



In most practical applications we deal with the multi-dimensional problems. For known 
solvable spherically symmetric potentials, the conventional WKB method does not re- 
produce the exact energy levels unless one supplements it with Langer-like correction 
terms. 

To overcome this problem in particular case of the Coulomb potential, a special tech- 
niques has been developed. In order for the WKBi approximation to give the exact eigen- 
values, the quantity + in the radial equation must be replaced by (/ + l/2)^ [|^. The 



reason for this modification (for the special case of the Coulomb potential) was pointed 
out by Langer (1937) [|l8l from the Langer transformation r = e^, U{r) = e^^'^X(x). 
However, for other spherically symmetric potentials, in order to obtain the appropriate 
Langer-like correction terms, another special transformation of the wf and its arguments 
is required. 

There are several other related problems in the semiclassical consideration of the radial 
Schrodinger equation, (i) The WKB solution of the radial equation is irregular at r — > 0, 

i.e. R^^^{r) oc r^/^/r, X = ^Jl{l + 1), whereas the exact solution in this limit is R{r) oc 
rK (ii) The radial equation has no the centrifugal term when / = 0, i.e. the radial problem 
has only one turning point and one can not use the WKB quantization condition derived 
for two-turning-point problems. However, solving the equation for Z > by known exact 
methods one obtains energy eigenvalues for all /. (iii) The WKB solution of the angular 
equation has analogous to the radial one, incorrect behavior at ^ — >■ 0: Qi^^^{9) oc 9^, 
Ij? = m? — h'^/A, while the exact regular solution in this limit is Q^{9) oc ^l™"'. Angular 
eigenfunction YQQ{9,ip) = const, i.e. no nontrivial solution exists. 

Practical use shows the standard leading-order WKB approximation always reproduces 
the exact spectrum for the solvable spherically symmetric potentials V{r) if the centrifugal 
term in the radial Schrodinger's equation has the form (/ + l/2)^^^/r^. As was shown in 
Ref. |ff9l the centrifugal term of such a form can be obtained from the WKB solution of 



the angular Schrodinger's equation if separation of the original three-dimensional equation 
has performed with the help of the correspondence principle. 

In Refs. IP, |l9l we have fulfilled the quasiclassical analysis of the three-dimensional 
Schrodinger's equation and suggested the quasiclassical approach for mult i- dimensional 
problems. In this approach the original wave equation is reduced to the form of the 
classical Hamilton- Jacobi equation without first derivatives. It was shown that, in the 
quasiclassical region the wave equation can be written in the form 



IS 



+ 



1 



+ 



4>{r} = 2m [E -V{r)]4>{r}. 



(32) 



Separation of this equation is performed with the help of the correspondence principle 
between classical and quantum-mechanical quantities. As a result of the separation, we 
have obtained a system of reduced second-order differential equations. 



-in4-] R{r) 
dr / 



M2 

2m {E-V{r)) 



R{r), 



(33) 
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■'4) 



d 



Ml 



e{9), 



(34) 



(35) 



where M^, M| are the constants of separation and, at the same time, constants of motion. 

Each of Eqs. (^3|)-(^) has the general WKBq solution of the form (||). Application 
of the quantization condition ([1^) to the angular equation (Q) results in the squared 
angular momentum eigenvalues P, [19| , 



(36) 



The same quantization condition appropriate to the radial equation (^) yields the 
exact energy spectra for all known solvable potentials and "insoluble" potentials, as well. 



7 Some practical applications 

To illustrate efficiency of the WKBq method consider several classic problems. The gen- 
eral solution (^ is the same for all problems in quantum mechanics. In case of multi- 
dimensional problems, we use the corresponding mult i- dimensional classical action. The 
oscillating part of the WKBq wf (|23| ) gives the asymptote of the exact solution of the 
Schrodinger's equation. Most complicate computational part in this approach is calcu- 
lation of the phase-space integral. In case of multi- dimensional problems, separation 
should be performed with the help of the correspondence principle between classical and 
quantum-mechanical quantities 0, |1^ . 

7.1 The linear harmonic oscillator V{x) = ^muj'^x^ 

Consider first the linear harmonic oscillator. Many problems in physics can be reduced 
to a harmonic oscillator with appropriate conditions. The general differential equation 
for oscillator potential can be solved using a technique that is frequently exploited in 
solving quantum mechanics problems. The eigenfunctions that are the solutions of the 
Schrodinger equation are the Hermite polynomials. In our approach, we directly ob- 
tain the asymptote of the exact solution and the corresponding well known exact energy 
eigenvalues. 

The quantization condition (|T9|) appropriate to the Schrodinger equation with the 
linear oscillator potential is 

J = y \j2mE — (muxydx = vr^ + . (37) 

The phase-space integral ( P7[ ) is calculated in closed form, / = E/lj, that results in the 
exact energy eigenvalues, E = ujh{n+l/2). The WKBq wf of the oscillator (the asymptote 
of the problem) is given by Eq. (^) and has the form of the standing wave. 
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IpO^x) = Cn COS ( ^PnX + ■ (3^ 



The normalization constant is calculated according to Eq. ( p^ with P„ = -\/2m£'„. 
7.2 The Coulomb problem V{r) = —a/r 

The Coulomb potential is a classic example of the exactly solvable problems in quantum 
mechanics. For this 2TP problem, the standard WKB method does not reproduce the 
exact energy levels unless one supplements it with Langer-like correction terms. In our 
approach, separation of the original 3D Schrodinger's equation with the help of the corre- 
spondence principle [|l9l results in the centrifugal term M^/r^ with given by Eq. (|36|) , 



which does not require any Langer type corrections. The quantization condition (|19]), for 
the Coulomb problem, is 



^ 2ma M2 , ^ / 1\ 

2mE + -dr = 27r/i + - , (39) 

r \ 2/ 

where the integral is taken about a contour C inclosing the turning points ri and r2. 
To calculate the integral (|39|) we use the method of stereographic projection. To do 
this, we should exclude the singularities outside the contour C, i.e. at r = and cxo. 



Excluding these infinities we have, for the integral (|39D, I = Iq + I^o, where Jq = —2'kM 



and loo = 2TTiam/ y/2mE. The sequential simple calculations result in the exact energy 
spectrum 



Note, for the energy of zeroth oscillations we have, from Eq. (|3), Eq = —a'^m/[2{h/2 + 
Mq)^], that apparently shows the contribution of the quantum fluctuations of the angular 
momentum, Mq = h/2, into the energy of the ground state Eq 0. The radial WKBq 
eigenfunctions, R^\r), inside the classical region [ri,r2] are written according to Eq. 
( p3[) and give the asymptote of the problem, 

R(^\r) = C„cos (^ip^r + , (41) 

The normalization constant, Cn, is calculated with the help of Eq. (P^, where P„ = 
2m\En\ and En is given by Eq. (^Of) . 

7.3 The isotropic oscillator V{r) = ^muj^r'^ 

The three-dimensional harmonic oscillator is another classic example of the exactly solv- 
able problems in quantum mechanics. As in case of the Coulomb potential, the WKB 
method does not reproduce the exact energy levels unless one supplements it with Langer- 
like correction terms. This 4TP problem is usually solved with the help of the replacement 
z = r^, which reduces the problem to the 2TP one. 
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This problem can be solved as the 4TP problem in the complex plane. Because 
of importance of the oscillator potential in many applications and with the purpose of 
further development of the WKB method, in Ref. [0 the problem has been solved by 
two methods, on the real axis as 2TP problem and then in the complex plane as 4TP 
problem. 

Consider the WKBq solution of the problem. In the complex plane, the problem has 
two cuts, between turning points ri, r2 and rs, r^. To apply residue theory for the phase 
space integral we need to take into account all zeroes of the wf in the complex plane, i.e. 
the contour C has to include both cuts. The quantization condition (|3T|) has the form 



I = f \^ 2mE - (mwr)2 - ^dr = 2^71 + ^) , (42) 

where /x = 4 and C is the contour about the cuts at r < and r > 0, respectively. The 
number = nr<o + ?^r>o, where nr<o and nr>o are the numbers of zeroes of the wf at 
r < and r > 0, respectively. For the 3D harmonic oscillator, because of symmetricity of 
the potential, we have nr<o = '^r>o = f^r, i-e. the total number of zeroes is = 2nr- 

Integral (|4^) is reduced to the above case of the Coulomb potential with the help of 
the replacement z = r"^. Integration result is I = tt{E/uj — M)/2 and we obtain, for the 
energy eigenvalues. 



En = UJ 



2h{nr + ^]+M 



(43) 



So far, as M = (/ + l/2)h, we get the exact energy spectrum for the isotropic oscillator. 
Energy of the ground state is £"0 = '^(^ + Mq), where Mq is the contribution of quantum 
fluctuations of the angular momentum M . The asymptotic eigenfunctions have the form 



(|38D, where the eigenmomentum P„ is calculated with the use of the energy eigenvalues 



m 



7.4 The Hulthen potential V{r) = -Vqc-''/''^ /{I - e 



The Hulthen potential is of a special interest in atomic and molecular physics. The po- 
tential is known as an "insoluble" by the standard WKB method potentials, unless one 
supplements it with Langer-like corrections. The radial problem for this potential is usu- 
ally considered at / = 0. However, in the approach under consideration, the quasiclassical 
method results in the nonzero centrifugal term at / = and allows to obtain the analytic 
result for all I. 

The quantization condition (Rl|) for the Hulthen potential is 



^2m^E + Vo^^^^j - -^dr = 2^h [n^ + - j . (44) 

In the region r > 0, this problem has two turning points ri and r2- The phase-space 
integral ( |4^ ) is calculated analogously to the above case. Introducing the new variable 



p = r/r-Q, we calculate the contour integral in the complex plane, where the contour 
C encloses the classical turning points pi and p2. Using the method of stereographic 
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projection, we should exclude the infinities at r = and oo outside the contour C. 
Excluding these infinities we have, for the integral (PI), / = /q + /qo, where Jq = — 27rM 
and /oo is calculated with the help of the replacement z = — 1 [H , 



2mrQ E 



Vn 



e-P 



e-P 



M2 , 

—dp 

p2 



-2nM + 2nroV^ 



m 



(45) 



Substituting the integration result into Eq. (|4^), we immediately get the exact energy 
spectrum 



En 



1 f2mVorl _ V 



imr, 



N 



(46) 



where = {rir + l/2)h + M is the principal quantum number. Setting in ( ^6]) M = 0, we 
arrive at the energy eigenvalues obtained from known exact solution of the Schrodinger's 
equation at / = 0. However, in our case M^m = Mq = h/2 aX I = and the principal 
quantum number is = [rir + l/2)h + Mq. As in the previous examples, this apparently 
shows the contribution of the quantum fiuctuations of the angular momentum into the 
energy of the ground state, Eq. 

7.5 The Morse potential V{r) = yo[e-2«(^Ao-i) _ 2e-«(V'^o-i)] 

The Morse potential is usually considered as one- dimensional problem at / = 0. In this 
case the problem has two turning points (note that the left turning point, ri, is negative) 
and can be solved exactly. In the general case, for Z > 0, we have an "insoluble" 4TP 
problem. 

For this potential, let us consider, first, the radial Schrodinger equation, which does 
not contain the centrifugal term at / = 0, 

(^-^^-^^ i?(r) = 2m[E-K)e"^"^"~''°^/''° + 2K)e~"(''""''^/''°]i?(r). (47) 
The quantization condition (|19|) appropriate to this equation is 

^ ^2m[E - •l/oe-2°('^-^o)Ao + 2Voe~''(''~''o)/ro]dr = 27Th (^rir + . (48) 

Introducing a variable reduce the phase-space integral to the well 

known one. Sequential simple calculations result in the exact energy eigenvalues 



En = -Vo 



1 - 



ah{nr + I' 



(49) 



roA/2mVo 

Now, let us deal with Eq. (^) for this potential, which contains the non-vanishing 
centrifugal term, h'^/Ar'^, at / = 0. In this case we have an "insoluble" 4TP problem. In 
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the complex plane, the problem has two cuts (z/ = 2), at r < and r > 0. Therefore, we 
apply the 4TP quantization condition ([31|), 



A2 



l2mr^[E - Voe-^''(p-^^ + 2V"oe-"(^-i)] - —dp = Anh 



rir 



(50) 



where we have introduced the new variable, p = t/tq. The contour C encloses the two 
cuts, but does not enclose the point r = 0. To calculate the integral (|^) we use the 
method of stereographic projection. For this, we should exclude the singularities outside 
the contour C, i.e. at r = and oo. Excluding these infinities we have, for the integral 

m i, 



27rro 



a 



V-2mE - y2mVoj , 



and for the energy eigenvalues this gives 



Er, 



a- 



2h{nr + i) + M 



H 2 



roA/2mT^ 



(51) 



(52) 



Setting in (|5^ ) / = 0, we obtain. 



En = -Vo 



1 



a[2h{nr 



Mn 



roA/2mV^ 



(53) 



Equation ( p3D for En is different from the expression (^) obtained from solution of radial 
Schrodinger equation for the Morse potential at / = 0. This difference is caused by the 
nonzero centrifugal term /Ar"^ in the radial equation (|33D at / = 0. Thus we obtain two 
results for the Morse potential: the known exact eigenvalues (^) obtained from solution 
of radial equation at / = and another result (p2D obtained from solution of Eq. (p3D for 
all /. 



7.6 The relativistic Cornell problem, V{r) = —a/r^nr 

In this paragraph, we use the WKBq approach to reproduce the exact energy spectrum for 
the famous funnel type potential (Cornell potential) [^| known as one which is "insoluble" 
exactly in terms of special functions. This potential is one of a special interest in high 
energy hadron physics, quarkonium physics, and quark potential models. Its parameters 
are directly related to basic physical quantities of hadrons: the universal Regge slope 
a' ~ 0.9(GeV/c)~^ of light flavours and one-gluon-exchange coupling strength for 



heavy quarkonia. A closer inspections reveals that all phenomenologically acceptable 
" QCD-inspired" potentials are only variations around the funnel potential. 

In relativistic theory, the Cornell potential represents the so-called "insoluble" 4TP 
problem. It is insoluble also from the viewpoint of the conventional 2TP WKB approxi- 



mation. Show that the WKBq method and the quantization condition (|3T|) derived above 



give the asymptote of the exact solution and yield the exact energy spectrum for the 
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potential. To solve the problem, we use the relativistic semiclassical wave equation ob- 
tained in Ref. [Q. For a two-particle system of equal masses, mi = 1712 = m, interacting 
by means of the Cornell potential, the relativistic radial semiclassical wave equation is 

{h = c = i) 



'-4 ^('■' 



/ a 

[m \- nr 

4 V r 



(54) 



where a = 40:^/3. The quantization condition (|31|) appropriate to Eq, 



IS 




a 



m 



+ Kr j 



47r n. 



(55) 



As in above examples, to calculate the phase-space integral (^) in the complex plane, 
we use the method of stereographic projection. Chose a contour C enclosing the cuts 
(and, therefore, zeros of the wf) at r > and r < between the turning points ri, r2 
and r3, r4, respectively. Outside the contour C, the problem has two singularities, i.e. at 
r = and 00. Using the method of stereographic projection, we have, for the integral 



5D, / = Jo + /oo, where 



-2'K^{l + l/2f + a^ 



(56) 



and the integral /oo is calculated in closed form with the help of the replacement of 
variable, z = 1/r, that gives 



27r 



a 



(57) 



Therefore, for E^, we obtain 



E' = 8 



K 



2{nr + 




(58) 



which is the exact result for the Cornell potential. 

It is an experimental fact that the dependence E^{1) is linear for light mesons (linear 
Regge trajectories [|^). However, at present, the best way to reproduce the experimental 
masses of particles is to rescale the entire spectrum assuming that the masses M„ of the 
mesons are expressed by the relation = E^ — C"^, where E^ = 8K{2nr + / + 3/2) is the 
exact energy eigenvalues for the linear potential and C is a constant energy (additional 
free parameter) ||2^. The equation for is used to shift the spectra and appears as a 
means to simulate the effects of unknown structure approximately. 

At the same time, the required shift of the spectra naturally follows from Eq. (|58|). The 
Coulomb- like term, —a/r of the Cornell potential takes into account the week coupling 
effect that allows to describe the Regge trajectories observed at experiment. For light 
mesons one may expect that the qq bound states will feel only the linear part of the 
potential, which gives the main contribution to the binding energy. We thus assume that 
the Coulomb term can be considered as a small perturbation. Then we obtain, from Eq. 

(I 
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El = 8k (^2nr + ^ - + ^) • (59) 

Equation (p9D does not require any additional free parameter and reproduces the linear 
Regge trajectories as it observed at experiment. We obtain the necessary shift with the 
help of the term —32Kas/3 which is the result of interference of the Coulomb and linear 
terms of the interquark potential. Note that we obtain the correct sign (minus) for this 
term only in the case of the scalar-like potential. Equation ( |55| ) describes the light meson 
trajectories with the accuracy ^1% pB, (see Table 1). 



Table 1. 


p family 


mesons 


State 




En'^P 




0.763 


0.768* 




1.319 


1.318* 




1.703 


1.691* 


13^4 


2.014 


2.037 




2.284 






1.703 


1.700* 




2.014 




2^Ds 


2.284 




2^F, 


2.525 





The symbol * denotes the masses of the fitted states; parameters of the fit are: = 
0.75 ± 0.03 and k = (0.14 ± 0.01) GeV^. 



8 Conclusion 

We have considered solution of the Schrodinger's equation in the zero order WKB approxi- 
mation. One might treat the WKBq solution as a very crude quasiclassical approximation. 
However, the WKBq solution represents the main term of the asymptotic series in theory 
of the second-order differential equations. A proper treatment of the WKBq approxima- 
tion opens the way to a simple asymptotic solution of the wave equation. 

It is a well known fact that the exact eigenvalues can be defined with the help of the 
asymptotic solution. The WKBq approach gives a simplest way to find the asymptote of 
the exact solution directly, without solution of the Schrodinger's equation, and results in 
quantization of the classical action. 

We have observed several important advantages of the zeroth approximation with 
reference to the commonly used first-order WKB approximation. Because the general 
WKBq solution (||) has no divergencies at the classical turning points, this allows a con- 
siderably simpler derivation of connection formulas and quantization condition. Using 
the general requirements of continuity and finiteness for WKBq function ipo{x) and its 
derivative ip'oix), we have derived a simple connection formulas, which allowed us to build 
the WKBq wf in the whole region and the corresponding quantization condition for the 
classical action. 
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If WKBi approximation is applicable at a distance from the turning points satisfying 
the condition |x — a;o| ^ A/47r, the WKBq solution considered in this work is applicable in 
the whole region. To demonstrate advantages of the WKBq approximation we have derived 
not only well known 2TP but, also, multi-turning-point quantization condition. This 
approach has been also used to derive the quantization conditions for periodic potentials 
and can be easily generalized for nonseparable problems. 

The most important advantage of the WKBq approach is that the WKBq quantization 
condition successfully reproduces the exact energy spectra for a// known solvable potentials 
and "insoluble" potentials, too. To demonstrate efficiency of the WKBq method, we have 
solved several classic problems. The method successfully reproduces the exact energy 
spectrum not only for solvable one-dimensional and spherically symmetric potentials but, 
also, for "insoluble" potentials with more than two turning points. 

For discrete spectrum, the WKBq eigenfunctions have written in terms of elementary 
functions in the form of a standing wave, AnCos{knX + The fact the eigenfunctions 
can be written in such form means the solution describes free motion of a particle-wave 
in the enclosure. This, in turn, means that the eigenmomentum is the constant of motion 
and the corresponding coordinate is cyclical. 

The zeroth approximation reproduces all basic results of quantum mechanics: eigen- 
values, asymptotes, probabilities, etc. Obviously, this is more than just an approximation. 
This is, at least, an alternative quantization method. What is more, we do not even need 
to solve the Schrodinger's equation to find the asymptote. After we have fixed the classi- 
cal action, using the existence/quantization formulas, we can write down the quantization 
condition for an arbitrary potential. 

This approach can be easily generalized for the non-separable problems. In this case 
the classical action can not be written as sum of actions for each variable (degree of 
freedom). As a result, we cannot have a unique quantum number for each degree of 
freedom. We will have one quantum number for the non-separable variables, and unique 
quantum number for each separable variable. Asymptote will have the same form, but 
the existence formulas (quantization condition) will be multi-dimensional integral over 
the non-separable variables. 

In the WKBq method, we use the same technique for all types of problems. The same 
simple rules formulated for 2TP problems work for /xTP problems, as well. In this sense, 
the WKBq approach can be considered as a general method for finding the asymptote of 
the Schrodinger's equation. 
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